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Abstract 

Non-local field theories which arise from p-adic string theories have 
vacuum soliton solutions. We find the soliton solutions at finite tem- 
perature. These solutions become important for the partition function 
when the temperature exceeds ms/g^ where is the string mass scale 
and go is the open string coupling. 



PACS numbers: ll.lO.Lm (Nonlinear or nonlocal theories and models), 11.27.+d 
(Extended classical solutions), 11.25.Sq (Nonperturbative techniques; string 
field theory), ll.lO.Wx (Finite-temperature field theory) 
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1 Introduction 



There is a special class of non-local field theories characterized by an infinite 
number of derivative terms in the form of an exponential. These field theories 
describe the open string tachyon in p-adic string theories [l]-[4]. In several 
recent papers we studied these non-local theories at finite temperature [5, 6]. 
These theories have no true particle degrees of freedom. All contributions to 
the equation of state arise from interactions. Perturbation theory can be used 
to study the thermodynamics up to a temperature on the order of ms/gl-i 
where itls is the string mass scale and Qo is the open string coupling. There- 
after the perturbative expansion breaks down and higher order terms become 
important. In addition, the vacuum energy density is positive and hierachi- 
cally suppressed with respect to the Planck scale, leading to the possibility 
that it may contribute to the cosmological constant. 

These non-local field theories are known to have vacuum soliton solu- 
tions. In this context soliton refers to a localized nonsingular solution to the 
classical field equation with finite action. They have an energy proportional 
to msl gl and therefore will contribute substantially to the partition function 
for temperatures of that order and higher. In this paper we study finite tem- 
perature classical solutions to the p-adic string theories. These solutions are 
non-analytic in the open string coupling Qo and therefore cannot be calcu- 
lated using perturbation theory. Hence they extend the results obtained in 
[5, 6]. 

The action for the p-adic theory is given by [l]-[4] 



S ^ 1}^ Id^'x 
9p 



I 

where □ = —dl -|- V|,_i in fiat space, and we have defined 



i0e-°/^V+^</>^+i 

r ^ p + v 



(1) 



^ and = . (2) 



9p glv-^ Inp 



The dimensionless scalar field 4>{x) describes the open string tachyon, is 
the string mass scale, defined by ml = l/2a' with a' the string tension, and 
go is the open string coupling constant. Though the action (1) was originally 
derived for p a prime number, in this paper we allow it to be any odd integer 
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equal to or greater than 3. For constant fields, the resulting potential takes 
the form 



U 



m 



D 



Its shape is shown in Figure 1. 
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Figure 1: (Color online) Inverted potential of the p-adic tachyon for p = 3, 
7 and p — t- oo. 

At finite temperature the scalar field must be a periodic function of imag- 
inary time r with period equal to the inverse temperature /3 = 1/T, namely, 
0(r + /3, x) = 0(r, x) [7]. The classical equation of motion is 



(4) 



where □ = d"^ jdr^ + V^. Substitution into the action gives S = —S with 



S 



1 f p — l\ m 



D 



2\p+l) gl Jv 



/3/2 
-/3/2 



rfr0P+^(r, x) 



(5) 



The quantity 5* is obviously non-negative for odd powers p. The classical 
equation of motion extremizes the action, which in this case means that it 
minimizes S. 
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All of the finite temperature calculations in our previous papers were 
done around the perturbative vacuum = 0, which is of course a trivial 
solution to the classical equations of motion. Generally, a nontrivial solution 
(ps will have an action with 5'^ > 0. If, for example, this solution is thought 
of as a soliton, and the centers of the solitons are widely separated in space- 
time, then they contribute to the partition function as a dilute gas. In the 
usual way [7] this leads to 

oo -| 

Zs^T.-i\(^VKexp{-Ss)Y (6) 

n=0 ^■ 

or 

\nZs^/3VKeM-Ss)- (7) 

Here K is the factor due to quantum fluctuations around the classical solu- 
tion. The factor I3V arises because the soliton could be centered anywhere 
within the spatial volume V and imaginary time interval /3. The nontrivial 
solutions are exponentially suppressed according to the actual value of their 
action relative to the trivial solution. Calculation of the factor K is beyond 
the scope of this paper. 

Since the equation of motion is nonlinear, there may be multiple non- 
trivial solutions at any given temperature, or none. Our goal is to find the 
solutions with the minimum value of S at finite temperature since the others 
will be exponentially suppressed in comparison. 

The outline of our paper is as follows. In Sect. 2 we derive differential and 
integral equations to be solved at finite temperature. In Sect. 3 we extend 
the diffusion equation method to finite temperature. In Sect. 4 we deduce 
from general considerations three types of solution at finite temperature. In 
Sects. 5 and 6 we calculate solutions at all temperatures which are even and 
odd in imaginary time, respectively. In Sect. 7 we calculate the value of 
the action for all these solutions to determine which are most important as 
a function of the temperature. In Sect. 8 we conclude with a summary of 
results obtained and potential future directions. 
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2 Differential and Integral Equations 

The classical equation of motion has soliton solutions, as first discovered in 
Euchdean space at zero temperature [3] . Recall that 

e-°/^V ^(ff. (8) 

This has the trivial solutions = 0, = 1, and = — 1 (if p is an odd 
integer). Now make the ansatz 

<t>{T,^)^h{T)h{x)h{y)h{z) (9) 

SO that each / satisfies 

It is easily verified that / must satisfy 

fix) = dx'e-^'^^'~-^'l'r{x') . (11) 

zyvr J-oo 

Once again, trivial solutions include f — 0, f — 1, and f — —1 (if p is an 
odd integer). A nontrivial solution is 

fix) = ±pi/2(f-i)e-(p-i)A^^-V4p . (12) 

Multiplying the /'s together results in a soliton in one, two, three, or four 
dimensions. 

For T > the scalar field must be periodic in r with period /3 [7]. The 
ansatz above is applicable, and the solutions in the three space directions are 
still valid. In the imaginary time direction the differential equation is still 

but now /(r + /9) = /(r). 

Any periodic function can be expanded in a Fourier series. 

oo 
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Here w„ = 2nTn is the Matsubara frequency. Substitution into the differen- 
tial equation results in 

oo 

/p(r)= £ c„e'^"/^'e^'^"^ (15) 

n=— oo 

On the other hand 

1^''' dTfnT)e-^"^^ - /3c^e'"-/^' (16) 

J-/3/2 

which implies that 

(3 1'2. 

c„ = e-"'/*^'r / drfPMe-'^-^ . (17) 

y-/3/2 

The integral equation at finite temperature is therefore 

oo f0/2 

f(r) = V Q-^i/M\ia:nrrp / dr'e-"^-^' (t') . (18) 

„_ „ ^-/3/2 



One may also expand /(r) in sines and cosines instead of complex expo- 
nentials. 



/(r) = Ao + ^ K cos(a;„r) + 5„ sin(a;„r)] (19) 

n=l 

In the usual way this leads to 

/.R 12 oo 
cir7^(T') + 2rVe-"'/^'x 
-/3/2 „=i 

cos(c^„t) / dT'P{T') cos(a;„r') sin(a;„r) / dr' P{t') sin(a;„T') 

7-/3/2 ^-/3/2 



(20) 

Q 

/e(T) = T X: e-"/^' cos(a;„T) / dr'/i'(r') cos(a;„T') , (21) 



This naturally separates into one integral equation if the solution is even 

r/3/2 

and another if the solution is odd 



oo 12 

Ut) = 2T ^ e--"/^' sin(a;.r) / dr' F,{t') sin(a;„T') . (22) 



n=l 



■/3/2 
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This separation also follows from eq. (18). 

The kernel of the integral equation (18) is just the thcta function of the 
third kind 



oo 



OMq)= E 9"e2""S (23) 



so that in more compact notation 

r/3/2 

-/3/2 

Now it is helpful to use the identity 



/(r) = T f^' dT% Unr - r'), e-(2'^^)^/^n /^(r') . (24) 



9, {u, e-^) = e-^y^'^e, e-V.^) (25) 



to write 



M = ^ E / dr'/^(rOe-^^(^'-^+^'^)^/^ . (26) 

This is an alternative to an expansion in terms of trigonometric functions. 
The expression (26) can be written as 

,M.JL £ /'-'^^.„,,(.„_„,,e--(.-.^^^ (27) 
where t„ = r' + ri/S, then as 

since f{r + /3) = /(/3), and finally as 

/(r) ^ ^ E /^""''^r7^(r')e--^(--)V^ 

^ ^ r dr7nr')e-^'(-'--)'/^ . (29) 

Thus we recover the general integral equation (11), but (26) is potentially 
more useful since it embeds the required periodicity of the solution explicitly 
while (11) does not. In the above derivation we assumed that the integral 
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and the summation commutes. This is justified since the infinite sum clearly 
converges, and the individual integrals are finite (as long as / is nonsingular) . 

In this section we have derived several integral equations which give solu- 
tions to the original differential equation. These equations may involve either 
trigonometric functions, such as eqs. (18), (20), (21) and (22), or Gaussians, 
such as eqs. (11) and (26). The original differential equation displays a sym- 
metry where / — )■ — )■ /, and changes sign. This symmetry may be 
applied to all of the integral equations too. 



3 Diffusion Equation 

Yet another approach is to cast the problem in terms of a diffusion equation 
[8, 9]. This can be accomplished by introducing an extra dimension labeled 
by the variable r. The original differential equation is written as 

exp(^-2r,^j/(T) = /nT), (30) 

where r* = 1/2M^. Now introduce a new function F(t, r) which satisfies the 
diffusion equation 

§,FM^^lF{r,r) (31) 

where 7 is an as yet unspecified constant which can be chosen later. This 
function must be periodic in r and can be expanded in trigonometric func- 
tions just as in eq. (19) except that now and depend upon r. Substi- 
tution into the diffusion equation determines those functions to be 

A-o{r) = Oo 

Ar,{r) = a„e-'^"''/^ n > 1 

Bnir) = 6„e-'^"'-/^ n>l (32) 

where the a„ and 6„ are constants. 

The relationship to the original ordinary differential equation is found by 
examining the operation 

exp(-2r*^j F{T,r) = 
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= E^-^(|:j F(r,r) = F(r,r-2r.7) (33) 

on account of the fact that F obeys the diffusion equation. Hence the op- 
eration of the exponetial operator in r is just a translation in the coordi- 
nate r. The desired solution / is related to F only at the point r*, namely 
/(r) = F(t, r*). Hence the proper equation of motion can be phrased as 

F(r,(l-27)r.) = Ff(r,r,). (34) 

Solutions to the original differential equation have been found by evolving 
some initial nontrivial configuration F(r, 0) to the final one F(r, r^,) via the 
diffusion equation. Although this is an interesting approach, so far we have 
not found any advantage over solution of the original differential or integral 
equations in one variable. 



4 General Considerations 

One can prove some general theorems analogous to the theorems for periodic 
solutions in real time [10]. Suppose we have a periodic solution satisfying 

/min < fir) < /max (35) 

with /mill < /max- Let US first assume that /max > 0. Then using eq. (11) we 
find that 

A/f roc 

fir) < ^ dr' e-^^(^'-)^/VLx = /^a. • (36) 

Taking the maximum value of /(r) on the left side we infer that /max < /max) 
which further implies that /max > 1- On the other hand, if /max < then we 
infer that — 1 < /max < 0. 

Now consider the lower limit of the oscillation. 

]\/f POG 

m >fr ^"'''^^'"^^'^V^in = /^in ■ (37) 

First assume that /min > 0. Taking the minimum value of /(r) on the left 
side we infer that /min > /mim which further implies that /min < 1- On the 
other hand, if /min < then we infer that /min < — 1- 
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Therefore, we conclude that there could exist three types of oscillations. 

(i) < /min < 1 and /max > 1- (u) /min < "1 and /max > 1- (hi) /min < "1 

and — 1 < /max < 0. These results are rather intuitive, and correspond to 
oscillations about = 1,0 = 0, and = — 1, respectively. Furthermore, note 
that there are conditions on the amplitudes of the oscillations, which are a 
result of the nonlinearity of the original differential equation. 

5 Even Solutions in Imaginary Time 

In this section we consider solutions that are even in r. First consider the low 
temperature limit when T <S M. Then the width of a vacuum soliton, 1/M, 
is much less than /3. This means that one can have a dilute gas of solitons. 
When — /3/2 < r < /3/2, only the n = term will contribute significantly in 
eq. (26). When r ~ /3, only the n = 1 term will contribute. When r ~ 2/3, 
only the n — 2 term will contribute. And so on and so forth. Then the 
periodic approximate solution to the integral and differential equations is 

oo 

J(^) ^ _t.pl/2(p-l) J2 e-{P-m''{r-pnf/ip _ ^3g^ 
n=— oo 

In the limit /3 — )■ oo with —13/2 <^ t <^ (5 /2 only the n = term survives 
and the vacuum sohition (12) is recovered. 

An alternative approach in the low temperature limit is to substitute the 
vacuum solution into eq. (21) and extend the limits of integration to infinity. 
This gives 

/(r) = ^^P"'"^'-"\P^,lf t e-^-^/(^-^)^^os(c.„r) 

= ±2pV2(p-i)yC^I,^3(7,rT,e-^(2.T)V(p-i)M^^ _ ^39^ 

Applying the identity (25) to this results in eq. (38). 

Since we have seen that oscillations can only exist around the points 
= 0, ±1, let us start by looking at small oscillations around = 1. To 
capture these let us make a truncated harmonic series expansion. 

/(r) =Ao + Ai cos{ut) (40) 



10 



Substituting (40) in the field equation (13), and using trigonometric identi- 
ties, we find that at the same level of truncation 

^p-l^ p(p-l) ^g-3^2 ^ ^ (4^) 

^^g-i^ p(p-l)(p- 2) ^,-3^, = e-VM^ (42) 
8 

In the lowest order approximation, -C 1, we find the usual harmonic 
oscillations around the minimum to be 

Ao = 1 

oj^ = M^lnp. (43) 

Because the solutions must be periodic in imaginary time with period /3 we 
must have cu = 2'KmT where m is an integer which, without loss of generality, 
we may assume to be positive. The smallest frequency corresponds to m = 1. 
This means that the solitons correspond to only a particular temperature 
Tc = M\/lnp/27r. How do we get solitons with other temperatures? This is 
where the anharmonicity comes into play and one realizes that the frequency 
changes as we change the amplitude of the oscillations. Going to the next 
order we have 

Ao = + (44) 

This means that the temperature decreases as the amphtude is increased. 

In the above analysis higher harmonics were neglected. The first overtone 
can be included by writing 

/(r) ^Ao + Ai cos(a;r) + A2 cos(2ujt) (46) 

Proceeding as before we find 

fir) = 1 - IpAl + A, cosiur) + -—^—-Aj cos{2u;t) + ■■■. (47) 

The expression for the frequency in terms of the amplitude is the same as 
before. 
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The solution that is even in r evolves with temperature in the following 
way. At T = the soliton solution is the Gaussian given by eq. (12). When 
T <S M a periodic solution is constructed from widely spaced Gaussians 
as expressed by eq. (38). As the temperature increases further the tails 
of the Gaussians begin to overlap and the solution becomes more uniform. 
Eventually it evolves into the form of cq. (46) with Aq ^ Ai ^ A2 and 
Aq ^ 1. The solution goes to the constant 1 at the well-defined critical 
temperature Tc — My/\np/27r; thereafter there is only the trivial solution 
f — 1. Exactly the same evolution happens for negative /. This behavior is 
in accord with the general considerations delineated earlier. 

Precise solutions for the full range of T can only be done numerically. 
It is convenient to use dimensionless variables t = T/M and u = Tr. The 
integral equation to solve, eq. (21), assumes the solution is an even function. 

fiu) = 2 / du'fPiu') + 4 V e-(2™*) cos(2™) / du' f^u') cos(2™') 
io io 

(48) 

There are no general methods for solving nonlinear differential or integral 
equations, and the fact that the differential equation under consideration 
has an infinite number of derivatives means that the problem is infinitely 
non-local. A straightforward iteration of the integral equation, starting with 
a low temperature or high temperature trial solution, does not converge. The 
reason is that the solutions always have a region where / > 1, and if the trial 
solution is too large by a small amount then raising it to a power p > 1 will 
take it even further away from the true solution. An important observation 
is to note a constraint that follows from eq. (48), which is 

/ duf{u) = / dur{u) . (49) 
JO Jo 

The method used here starts with a trial solution and iterates, where after 

each iteration, the new solution is scaled by a constant to satisfy this con- 
straint. In this way convergence is achieved quite rapidly. Typically we stop 
the iteration after an accuracy of at least eight significant figures is attained. 

Some illustrative results are plotted in Fig. 2 as a function of Tr, and 
in Fig. 3 as a function of Mr. Referring to Fig. 2, the profile is described 
very well by the Gaussian for T < 0.5Tc. It then transitions to a cosine 
before fiattening to the constant 1 for T > T^. This is exactly the behavior 
predicted by the analytical approximations. 



12 




Tt 



Figure 2: (Color online) Even soliton profile for p = 3 as a 
function of imaginary time r for different temperatures: T/Tc = 

0. 01, 0.10, 0.30, 0.50, 0.75, 0.85, 0.95. For T/T^ > 1 the solution is exactly 

1. As the temperature increases, the solution transitions from a narrow 
Gaussian, to a cosine wave, to the constant 1 at T^. 
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Figure 3: (Color online) Even soliton profile for p = 3 as a 
function of imaginary time r for different temperatures: T/Tc = 
0.10, 0.30, 0.75, 0.85, 0.95, 1.00, from top left to bottom right, in units of 
the p-adic string scale M. 
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6 Odd Solutions in Imaginary Time 



Now consider an oscillation that is odd in r. Such an oscillation can only 
occur about = 0. The harmonic expansion starts out as 



/(r) = Bi sin(a;T) + B3 sin(3a;T) . 



(50) 



In order that the solution have the required periodicity implies that a; = 
27r(2m + 1)T where, without loss of generality, m = 0, 1, 2, .... Substitution 
into the differential equation, and matching only the sin(a;T) term, results in 



Bi 



(p + 1)!! 
2p\\ 



i/{p-i) 



,a;V(p-l)M2 



(51) 



Matching both the fundamental and first overtone results in 



Bi 
B3 




1 ip+^r 



^a;2/(p-l)M2 



P ~ 1 ) „-8a;VM2 
.P+V 



(52) 



This solution has the characteristic feature predicted previously; namely, the 
amplitude of oscillation exceeds 1. In fact, the amplitude grows exponentially 
as T^/M^, and the fundamental frequency dominates at large temperature. 
Although solutions exist for any integer value of m, the solution with m = is 
the most important because the others have actions which are exponentially 
larger. 

When T > M/An the first term in the series (50) dominates. As T 
decreases the amplitude decreases, and more and more terms in the series 
become important. There is no known nontrivial solution at T = 0, so it is 
interesting to see what happens in the limit T ^ 0. In fact, one can construct 
a solution both periodic and odd in imaginary time by adding the vacuum 
Gaussian solutions in the following way. 



/(r)=pV2(P-i) J2 {-ire- 



(p-l)M2(T-/3(n+l/2))2/4p 



(53) 



When T < M/An, these Gaussians are spaced much further apart than their 
widths, and so the differential and integral equations are satisfied to very 
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good approximation. In particular, there is a positive Gaussian centered at 
T — (3/2 and a negative Gaussian centered at r = — /3/2. As T increases 
these naturally go over to Bi sin(27rT'T). As T goes to zero, the spacing in r 
between adjacent Gaussians diverges as (5 = 1/T . 

Once again, precise solutions at intermediate temperatures require nu- 
merical calculation. For an odd solution the equation to solve is 

°° (•1/2 

f{u) = 4 ^ e-(2""*)' sin(2™) / du' f{u') sin(2™') . (54) 

n=l 

In this case the constraint we use takes the form 

rl/2 1-1/2 

/ dii/(ii) sin(27rii) = e-^^'^*) / duf^{u)sin{2T:u) . (55) 
Jo Jo 

As with the even solutions, convergence is achieved rather rapidly. 

Some illustrative results are plotted in Fig. 4 as a function of Tr, and 
in Fig. 5 as a function of Mr. Referring to Fig. 4, the profile is described 
very well by a pair of Gaussians for T < O.ST^. It then rapidly transitions 
to a sine with increasing amplitude for T > Tc. This is exactly the behavior 
predicted by the analytical approximations. 



7 Action 



As discussed in the introduction, the contribution of solitons to In Z decreases 
exponentially with their action. The spatial distribution of the solitons is 
always given by the Gaussian of Eq. (12). Inserting this into the action gives 



S{T) 



1 



p 



2 \P+iJ gl 

where the dimensionless integral 



— 1 



(D-l)/2 



(56) 



/(T)=m, / dTr+\T) 



(57) 



determines the temperature dependence. 
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-0.5 -0.25 0.25 0.5 



Tt 

Figure 4: (Color online) Odd soliton profile for p = 3 as a function of imag- 
inary time r for different temperatures: T/Tc = 0.01, 0.10, 0.50, 1.00, 1.50. 
As the temperature increases, the solution transitions from a pair of narrow 
Gaussians to a sine wave whose amplitude increases with temperature. 
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Figure 5: (Color online) Odd soliton profile for p = 3 as a function of imag- 
inary time r for different temperatures: T/Tc = 0.10, 0.50, 1.00, 1.50, from 
top left to bottom right, in units of the p-adic string scale M. 
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First consider the low temperature limit. For T < 0.5Tc the even solution 

gives 



to good approximation. Meanwhile the odd solution gives 



7(T) = 2, ^p^/M 



p2 - 1 



(58) 



(59) 



to good approximation, which is twice the value of the even solution. This is 
because there are two Gaussian distributions inside the interval of integration. 
See Figs. 2 and 4. 

Now consider the high temperature hmit. For T > Tc the exact result 
for the even solution is 

m = ^ . (60) 



T 



For T > Tc the odd solution gives 



/(T) 



V2T 



(p + 1) 
2p\\ 



1 2/(p-l) 



exp 



1 j V M y 



(61) 



which increases exponentially in T'^/M'^. 

We can compute the contributions from the classical soliton configura- 
tions numerically. These are displayed in Fig. 1. As discussed above, the 
high temperature behaviors are opposite. The action for the even soliton 
solution decreases as 1/T for T > Tc while the action for the odd soliton 
solution grows exponentially. At low temperature the action for the odd so- 
lution is twice that of the even solution. From these computations, we can 
deduce that the odd soliton contribution is unimportant at all temperatures, 
whereas the even soliton contribution is important when T > nis/g^- 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 

T/T, 



Figure 6: (Color online) Contribution /(T) to the Euclidean action of the 
even (top) and odd (bottom) solitons for p = 3 as a function of T/Tc. The 
green dashed lines show the asymptotic behaviour at low and high tempera- 
tures. 
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8 Discussion 



In this paper we have found finite temperature sohton solutions to a class of 
non-local field theories arising from p-adic string theory. Analytic solutions 
were derived for both low and high temperatures, delineated by a critical 
temperature Tc ~ iJis/ g1- In the intermediate temperature region numerical 
solutions were readily found. 

These soliton solutions have finite action even at zero temperature. This 
means that they will contribute to the vacuum energy density or cosmo- 
logical constant. However, their contribution is suppressed exponentially as 
Q-cp/9o ^ where Cp is a number dependent upon p, relative to the perturbative 
contribution [6]. Hence they will be unimportant for open string couplings 
^o< 1- 

In order to quantify the soliton contributions to the partition function, 
quantum fluctuations around the classical solutions must be calculated. That 
work is in progress. 
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